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ABSTRACT 

We present Monte Carlo models of open stellar clusters with the purpose of mapping out 
the behavior of integrated colors with mass and age. Our cluster simulation package allows 
for stochastic variations in the stellar mass function to evaluate variations in integrated cluster 
properties. We find that UBVK colors from our simulations are consistent with simple stellar 
population (SSP) models, provided the cluster mass is large, Mduster > IO^Mq. Below this 
mass, our simulations show two significant eflFects. First, the mean value of the distribution of 
integrated colors moves away from the SSP predictions and is less red, in the first 10^ to 10^ 
years in UBV colors, and for all ages in {V — K). Second, the Icr dispersion of observed colors 
increases significantly with lower cluster mass. The former we attribute to the reduced number of 
red luminous stars in most of the lower mass clusters and the later we attribute to the increased 
stochastic effect of a few of these stars on lower mass clusters. This later point was always 
assumed to occur, but we now provide the first public code able to quantify this effect. We are 
completing a more extensive database of magnitudes and colors as a function of stellar cluster 
age and mass that will allow the determination of the correlation coefficients among different 
bands, and improve estimates of cluster age and mass from integrated photometry. 
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Introduction 



When dealing with the integrated light from ex- 
tremely massive systems, such as entire galaxies, 
one can reasonably assume there are enough stars 
in the system to represent all critical stages of evo- 
lution, regardless of stellar mass or age. When 
studying very massive stellar clusters, 100,000 Mq 
or more, while all members are (typically) of the 
same age, the isochrones will be fully p opulated at 
all ages (e.g. iBruzual fc Charlot|[2003l) . 

However, typical open star clusters do not have 
such high masses. Two open star clusters of nearly 
identical age can not be expected to have the 



open clusters and associations: 



same integrated light properties. This is because 
of natural cluster-to-cluster variations in the ex- 
act stellar masses found within the cluster, which 
give way to the cluster's integrated ligh t proper- 
ties ( e .g. Cervino fc Valls-Gabaud 20091 : Buzzoni 
19891: IChiosil Il989l l In their final evolutionary 



states, stars become very luminous and/or possess 
rather extreme red or blue colors. This can cause 
observable deviations in integrated colors of the 
cluster as compared to the integrated colors com- 
puted in the high mass limit. This effect is partic- 
ularly great in young to middle-aged open clusters 
where their dying stars are the most luminous and 
have the most extreme and rapid color changes. 
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This poses a serious problem to astronomers be- 
cause integrated colors are the primary way in 
which the age and m ass of unresol ved stellar clus- 
ters are derived (e.g. Larsen 2009() . 

The fact that integrated colors are affected by 
the finite value of cluster mass and the fluctua- 
tions in the IMF ha s been known and studied by 
many authors (e.g. lBruzual|l2002t lBruz ual"2009l: 
Ceryino &: LuridianaT 2004 Cerviho &: Luridiana 



explore the colors and color ranges of a stellar clus- 
ter as the mass drops below this limit, introducing 
mass as a variable in this investigation. 
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2000; La ncon fc Fo uesneau 2009; Foues neau fc Lancon 
2009; Brocato et al... .1999: Brocato et al.l l2000t 



Buzzonilll989l: Chiosil[l989[ ICantielb et al'l 



2003; 



Raimondo et al.'2005; Raimondo 2009; Piskunov et al 
2009; Santos & Fr ogel 1997; Gonzalez et al., 2004t 



Gonzalez-Lopezli ra et al.ll2005t iGonzalez-Lopezlira et al 
20101 : iMafz ApeMnid 120091 : IPessev et al.l l2009l ). 



We have computed the mean value of {B — V)o, 
{U — B)o and {V — K)o integrated colors as a 
function of mass and age for a simple stellar clus- 
ter. Our results - MASSCLEAN co/ors database - 
are based on over 30 million Monte Carlo simula- 
tionH. Broadly speaking, the simulations provide 
UBVRIJHK magnitudes as a function of age 
and mass. For this demonstration we present just 
UBVK, for two metalhcities, Z = 0.019 (solar) 
and Z = 0.008 (Large Magellanic Cloud, LMC) 
and employ the two most commonly used st ellar 
evolutio n models, Padova ('Marigo et al.l[2008 ^ and 
Geneva (jLejeune fc Schaerc r 2001,). The simula - 



2009t ). This is a pubhc and open source code. 



We present a large number of Monte Carlo simu- 
lations, which constitute the MASSCLEAN colors 

database, using the MASSCLEAN (MASSiv e CLuster . , . .^^ „„„ . 

„ -, ^ . J .„ -, • \ 1 rn,Ti TTr 1 tions were done usmg Kroupa IMl* ( Kroupa 2002 ) 

Evolution and ANalvsis) packagd Popescu &: Hanson . 7^ 

^ ^ — ^ ^ with 0.1 M© and 120 Mq mass limits. The age 

range of [6.6,9.5] in log{age/yr) was choosen to 
accomodate both Padova and Geneva models. An 
average of 5000 clusters were simulated for each 
mass and age. These results are presented in Fig- 
ures 1 through 5. Comparison with observational 
data from LMC and Milky Way clusters is dis- 
cussed in [231 



While it was first designed to simulate particular 
clusters, new additions to the package allow it to 
run in Monte Carlo mode. Also, all of the tools 
used to perform statistical analysis described in 
this work are included, so we are able to directly 
predict the mean value of the distribution of colors 
and color dispersion as a function of cluster mass 
and age. 

MASSCLEAN is a new stellar cluster image and 
photometry simulation package. It uses an algo- 
rithm which populates the simulated cluster with 
a discrete number of tenable stars and then evolves 
e a,ch individual star f ollowing isoch rones (Padova 
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The "blueing" of mean integrated col- 
ors for lower mass clusters 



Marigo et al.l l2008l or Geneva - ISchaller et al 



I992I : Mevnet et al. 1994 with isochrones obtained 



using Basel a tmosphere libraries as presented in 



Lejeune fc Schacrer 200 ll) . Integrated photometry 



of the simulated cluster is derived by simply sum- 
ming up the stellar flux at any point in the cluster 
evolution. We have tested that at very high cluster 
mass, the "infinite mass limit" of Mduster > 10^ 
Mq , the integrated colors are consistent with stan- 
dard simple stellar population (SSP) models when 
using identical clust er properties of age, IMF, evo- 
lutionary code, etc. ( Popescu fc Hanson '2009). In 
this Letter, we present results using MASSCLEAN to 



Figures ID through m are all similarly presented. 
The mean value for {B — V)o, [U — i?)o and 
{V — K)q distribution of integrated colors as gen- 
erated by MASSCLEAN are given in the top of each of 
the three panels in each figure. In Figures [T] and 
[21 we use the Padova evolutionary models with 
log{age/yr) step of 0.01, solar and LMC metallic- 
ity, respectively. In Figures [H and IH we present 
the mean values using the Geneva evolutionary 
models with log(age / yr) step of 0.05, solar and 
LMC metallicity, respectively. The difference in 
smoothness seen in the Padova and Geneva sets is 
due to a smaller log{age/yr) step in the Padova 
case (0.01, versus 0.05 for Geneva). 



^http : //www . physics . uc . edu/- popescu/massclean/ | 



■^corresponding to over 200, 000 MASSCLE AM runs (for a com- 
plete description of a MASSCLEAN run see lPopescu fc HansonI 
l2009h 
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Fig. 1. — MASSCLEAN integrated colors for Padova models Z — 0.19 using 0.01 log{age/yr) step. The thick line 
corresponds to the infinite mass limit (10^ Mq in our simulations), (a) Upper panel: {B — V)o as a function of mass 
and age; (b) Middle panel: {U — B)o as a function of mass and age: (c) Lower panel: {V — K)o as a function of mass 
and age. The Milky Way clusters are from iHancock et al.l (|2008l ) and iKharchenko et al.l (|2009l ). 
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Fig. 2. — MASSCLEAN integrated colors for Padova models Z — 0.008 using 0.01 log{age/yr) step. The panels are the 
same as in Fig. 1, except for the data points. Here, the clusters are from the LMC and found in iHunter et al] l|2003l ) 
and Marigo et al. (2008). 



4 



0.75 




6.75 7.00 7.25 7.50 



7.75 8.00 8.25 
log(age/yr) 



8.50 8.75 9.00 9.25 



9.50 



1o 100000 Ms„„ 
1o 50000 Ms„„ 
1c5 30000 Ms„_ 
1o 15000 Ms„. 
1o 10000 Ms„. 
la 7000 Ms,i„ 



iobOOO Ms„„ 
la 3000 Ms„- 
1a 2000 Ms„. 
1a 1500 Ms„. 
la 1000 Mgun 
1a 700 Mq,,„ 



la 500 Mgun 
Geneva Z=.020 
Hancock et al. (2008) 
Kharchenko et al. (2009) 



(a) (B - V)o 




6.75 7.00 7.25 7.50 



7.75 8.00 8.25 

log(age/yr) 
(b) {U - B)o 



8.50 8.75 9.00 9.25 



9.50 



3.50 
3.00 
^ 2.50 
^ 2.00 
— 1.50 
^ 1.00 

> 0-50 
0.00 
-0.50 

4.00 

„ 3.00 

I 2.00 

^ 1.00 

> 00 
-1.00 
-2.00 




Geneva Z=.020 - 
Geneva Z=.02O, 500-1 00000 Ms„. 
i i I i i i i I i i i i I i i i i I 




Geneva Z=.020 
Kharchenko et al. (2009) 



6.75 7.00 7.25 7.50 



7.75 8.00 8.25 

log(age/yr) 
(c) {V - K)o 



8.50 8.75 9.00 9.25 



9.50 



Fig. 3. — MASSCLEAN integrated colors for Geneva models Z = 0.20 using 0.05 log{age/yr) step. The panels are the 
same as in Fig. 1, including the Milky Way clusters shown. 
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Fig. 4. — MASSCLEAN integrated colors for Geneva models Z = 0.008 using 0.05 log{age/yr) step. The panels are the 
same as in Fig. 1, except for the data points which are from the LMC. 
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The upper graph in each of the 3 panels (a) of 
Fig. [T] through 2] show the mean value of the inte- 
grated cluster color as a function of age. Thirteen 
masses are presented starting at 10^ Mq (dark 
solid line, and virtually the same as the SSP codes 
predictions) and going down to 500 M© (dotted 
lines). The dotted lines show a trend away from 
the SSP predicted values to bluer mean colors, as 
the cluster mass decreases for young clusters up to 
an age of as much as 10* years. This occurs with 
both solar and LMC metallicities and both evolu- 
tionary models, Geneva and Padova. More strik- 
ing, as the mass decreases our simulations show 
a blueing of the integrated {V — K)o color for all 
cluster ages simulated (bottom panel). 

2.2. The increased range in the observed 
integrated colors of low mass clusters 

The color dispersion (presented as Icr - stan- 
dard deviation - about the mean) for different val- 
ues of mass in the 500 — 100, 000 Mq interval is 
presented in the lower graphs of each panel. Here 
we use differing levels of grayscale tone to repre- 
sent the cluster masses. The most massive clus- 
ters show their la range in the darkest grayscale, 
while lower mass clusters are shown with increas- 
ingly lighter shades. Recall, the mean value of the 
distribution of integrated colors moves to bluer col- 
ors as the cluster mass decreases. This explains in 
part the greater color range seen in the low mass 
clusters on the blue side. 

The expected number of stars in any given 
evolutionary phase is fractional. This fact is not 
relevant if the cluster has enough stars (is massive 
enough) so variations oi +/— one in a given evo- 
lutionary phase have no impact. However, when 
the mass of the cluster decreases, this number 
becomes lower than one in the red, luminous evo- 
lutionary phases, so most of the clusters in the 
simulation will have no luminous red stars at all, 
since the number of stars in a given phase must 
be an integer. This produces a bluer mean color. 
However, a few clusters will have one (or a few) 
luminous red stars, so the result will be an excess 



2004 : ICerviho fc Luridianal 120061 : iFagiolini et al 
20071 ). While the correlation between the mean 



in bo t h sides (blue and red) (ILancon & Mouhcine 
2002 ; ICerviho fc Luridianal l2004t IFagiolini et al. 
20071 ). In this situation it is difficult to ob- 
tain a clear correlation of a particular (real) 
clust er integrated color w" i th mass for a single 
age ( Santos &: Frogel 1997 : Cerviho &: Luridianal 



value (of the distribution of integrated colors) and 
mass is apparent in our plots, an accurate quanti- 
tative description could be obtained only using a 
larger number of simulations. However, we hope 
the predicted increased color range with lower 
mass should be a strong enough signal with a re- 
alistically observable number of clusters. In order 
to verify this results with observations, perhaps 
100s of clusters would be needed and they must 
have well constrained mass and age. 

2.3. The effect of low cluster mass on the 
UBV color-color diagrams 

The color-age diagrams of Fig. |T] through H] 
are not typically used by researchers because age 
is often not known for unresolved stellar clus- 
ters. In fact, the great utility of SSP models 
has been to provide a way to derive cluster ages 
based on the easily observed integrated magni- 
tudes and their location in a color-col or diagram 
(ISearle et all ll980llGirardi et al.l Il995l) . In Fig- 
ure [5] we provide such diagrams. Here we have 
re-plotted the simulations already shown in Fig. 
[T] through |4l including the mean color lines and 
the greyscale la ranges of those colors, just like 
our earlier figures. Here the four panels repre- 
sent the two evolutionary models (Padova and 
Geneva), with each of those using LMC and Milky 
Way metallicities. We include the same real data 
shown in Figures |T] through 2] for the LMC and 
Milky Way. Again we provide both mean colors 
as a function of age, and the la range of that 
mean color. A similar color scheme for the range 
is employed, the darkest grayscale representing 
the most massive cluster simulated, while lighter 
shades represent lower mass clusters. It is easy 
to see the very large dispersion of the lower mass 
clusters, particularly in the blue, young portion of 
Fig. [5l While its rather crowded, we've attempted 
to demonstrate the location of the mean color as 
a function of mass through the use of differing line 
styles in the plot. As one might expect from our 
four earlier figures, most of the predicted mean 
variation occurs in the upper left portion of the 
diagram, corresponding to the youngest clusters. 
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Fig. 5. — MASSCLEAN integrated UBV color-color diagrams, (a) Using Padova models Z = 0.19 and Z = 0.008 and 
(b) using Geneva models, Z — 0.20 and Z = 0.008. The thick line corresponds to the infinite mass limit (10^ Mq 
in our simulations). On top is given Milky Way clusters from lHancock et al.l (|2008l ) and the bottom panels provide 
data for LMC clusters, taken from iHunter et aP l|2003l) . 
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2.4. Comparison with real cluster data 

In Fig. [Uth r ouKh H we p rovide data tak e n from 



Marigo et al 



(I2008D and 'Hunter et al.' (200^, 
for the LMC and iKharche nko ct al. (2009), and 
iHancock et"al1 (|2008l ) fo r the Milky Way. In Fig. 
[51 we provide data from Hunter et al.l ( 2003f ). for 
the LMC and lHancock et al.l (|2008[) for the Milky 
Way only. The other datasets did not include 
[/-band. In Fig. [SJ there appears to be as large 
a variation of observed colors along the portion 
of the color-color diagram for the LMC clusters 
where old clusters are expected to be found (in 
the lower right) as young clusters (in the upper 
left). This seems incongruous with the rather tight 
fit of panels (a) and (b) in Figures [2] and [H Also, 
in Fig. [5l we see that the youngest clusters from 
the LMC are found to lie almost entirely within 
the gray scale error range when they do not do so 
in Fig. [T] through 131 The first is due to a selection 
effect of the data. Clusters deviating greatly from 
the expected SSP curve in the UBV color-color 
diagram of Fig. |5l are not typically aged and 
thus will not appear in Fig. |2l and IH cleaning up 
those diagrams considerably. To understand why 
the younger clusters lie almost entirely within the 
grayscale of Fig. [5] and not Fig. [1] through |4l one 
needs to recognize that in the blue range, the lines 
representing the mean double back on themselves, 
making the region very complex. A quick check 
shows that the greatest range of grayscale in Fig. 
[5l never extends beyond the greatest range seen in 
Fig. [H through m 

Looking at the figures one notes the data does 
not appear to be distributed in a manner consis- 
tent with our la ranges. In fact, one could say, 
the colors of real clusters shown often extend to 
values much redder than our la range given in 
the figures (more often than la would indicate for 
the data). This is because of the way in which we 
have shown the color range with mass. To provide 
a la range is most useful when the distribution 
is Gaussian. Unfortunately, the distribution of 
observed colors is far from a well behaved Gauss 
about the mean for early ages and low mass clus- 
ters. The distribution of colors with age in Fig. 
[I] through Ul becoming increasingly bimodal as 
mass decreases. Thi s is consistent with the fluc - 
tuations d escribed by'Lancon & Mouhcine' 



Cerviho &: Valls-Gabau d (2003); Ccrviho & Luri diana 
(|2004l ): ICervifio fc Luridianal (|2006t ) . Wewihshow 



that our simulations are consistent with these very 
red colors observed in low mass systems (Popescu 
& Hanson 2010, in prep.). We are running ad- 
ditional simulations and eventually will present 
more than 50 million clusters to better constrain 
the exact behavior of the colors as a function of 
mass and age. 

3. Discussion and Conclusions 

While it is entirely expected that integrated col- 
ors computed by the modern SSP codes in the in- 
finite mass limit will work best when applied to 
very massive systems, its been difficult to ascer- 
tain the error involved when one applies their pre- 
dictions to lower mass systems. Our new stellar 
cluster simulation routine, MASSCLEAN, provides a 
means to determine this through Monte Carlo sim- 
ulations of realistically modeled stellar clusters. 

We used MASSCLEAN to simulate 30 million 
stellar clusters over a mass range typical of 
open stellar clusters and for ages in the [6.6,9.5] 
log{age/yr) range. Our results indicate that as 
cluster mass decreases, the mean value of the dis- 
tribution of colors in the first tens (to hundreds) 
of millions of years will not reach to such high red 
values as predicted by SSP models. Most of this 
very red color is brought on by very luminous red 
stars. Lower mass clusters will have fewer highly 
luminous red stars to pull the colors so red. How- 
ever, the stochastic variation of stars in these very 
luminous, quickly evolving phases, and thus very 
red colors, makes it still possible just not so likely. 
Thus we also find the typical range of colors ob- 
served increases significantly as the cluster mass 
decreases. This is because the color depends so 
strongly on the most luminous stars and a single 
star has the potential to produce greater devia- 
tions in the color of smaller clusters. Our results 
indicate a bluer mean value of the distribution of 
colors for lower mass clusters compared to SSP 
predictions in the infinite mass limit. 

We continue to simulate still more clusters and 
expand the MASSCLEAN co/ors database to pro- 
vide a complete sampling of colors and absolute 
magnitudes as a function of cluster age and mass 
(the work presented in this Letter was done dur- 
ing a six month period, using multiple comput- 
ers). When complete, this database will pro- 
vide a probabilistic grid from which a statistical 
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inference code will be based to aid researchers 
in deriving cluster characteristics from their ob- 
served integrated photometry. We have completed 
some preliminary work to provide proof of concept 
( Popescu fc Hanson 2010l ). With additional sim- 
ulations (Popescu & Hanson 2010, in prep.), we 
will be able to obtain the correlation coefhcients 
among different bands. 

We are grateful to suggestions made to an early 
draft of this work by Rupali Chandar and Deidre 
Hunter. We thank the referee for useful comments 
and suggestions. This material is based upon 
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